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ABSTRACT 

<N . 

^ ■ We examine the relation between the galaxy cluster mass M and Sunyaev-Zeldovich (SZ) effect signal D 2 A 7 50 o for a sample of 19 objects for 
which weak lensing (WL) mass measurements obtained from Subaru Telescope data are available in the literature. Hydrostatic X-ray masses are 
\l derived from XMM-Newton archive data, and the SZ effect signal is measured from Planck all- sky survey data. We find an M WL -D 2 A Y 500 relation 
that is consistent in slope and normalisation with previous determinations using weak lensing masses; however, there is a normalisation offset 
vn ■ with respect to previous measures based on hydrostatic X-ray mass-proxy relations. We verify that our SZ effect measurements are in excellent 
-==j= agreement with previous determinations from Planck data. For the present sample, the hydrostatic X-ray masses at R 50 o are on average ~ 20 per 
^ cent larger than the corresponding weak lensing masses, which is contrary to expectations. We show that the mass discrepancy is driven by a 
■ difference in mass concentration as measured by the two methods and, for the present sample, that the mass discrepancy and difference in mass 
f== ^l concentration are especially large for disturbed systems. The mass discrepancy is also linked to the offset in centres used by the X-ray and weak 
• • 1 lensing analyses, which again is most important in disturbed systems. We outline several approaches that are needed to help achieve convergence 
^ , in cluster mass measurement with X-ray and weak lensing observations. 

, Key words, cosmology: observations - galaxies: clusters: general - galaxies: clusters: intracluster medium - cosmic background radiation - 
?-H ■ X-rays: galaxies: clusters 
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1. Introduction 

Although the Sunyaev-Zeldovich (SZ) effect was discovered in 
1972, it has taken almost until the present day for its potential to 
be fully realised. Our observational and theoretical understand- 
ing of galaxy clusters has improved immeasurably in the last 
40 years, of course. But recent advances in detection sensitiv- 
ity, together with the advent of large-area survey capability, have 
revolutionised the SZ field, allowing vast improvements in sen- 
sitivity and dynamic range to be obtained (e.g., Pointecouteau 
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et al. 1999; Komatsu et al. 1999; Korngut et al. 2011) and cat- 
alogues of tens to hundreds of SZ-detected clusters to be com- 
piled (e.g., Vanderlinde et al. 2010; Marriage et al. 201 1 ; Planck 
Collaboration VIII 2011; Reichardt et al. 2012). 

The SZ signal is of singular interest because it is not affected 
by cosmological dimming and because the total SZ flux or in- 
tegrated Compton parameter, 7sz, is expected to correlate par- 
ticularly tightly with mass (e.g., Barbosa et al. 1996; da Silva 
et al. 2004; Motl et al. 2005; Nagai 2006; Wik et al. 2008; 
Aghanim et al. 2009). SZ-detected cluster samples are thus ex- 
pected to range to high redshift and be as near as possible to 
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mass- selected, making them potentially very powerful cosmo- 
logical probes. Notwithstanding, a well-calibrated relationship 
between the total mass and the observed SZ signal is needed to 
leverage the statistical potential of these new cluster samples. 

In fact the relationship between mass and Ysz is still poorly 
determined, owing in large part to the difficulty of making suf- 
ficiently precise measurements of either quantity. Moreover, the 
majority of mass measurements used to date (e.g., Benson et al. 
2004; Bonamente et al. 2008; Andersson et al. 2011; Planck 
Collaboration XI 2011) have relied on X-ray observations that 
assume hydrostatic equilibrium, which many theoretical studies 
tell us is likely to result in a mass that is systematically underes- 
timated by about 10-15 percent due to neglect of bulk motions 
in the intracluster medium (ICM; e.g. Nagai et al. 2007; Piffaretti 
& Valdarnini 2008; Meneghetti et al. 2010). This effect is now 
commonly referred to in the literature as the "hydrostatic mass 
bias". 

In this context, weak lensing observations offer an alterna- 
tive way of measuring the total mass. As the weak lensing ef- 
fect is due directly to the gravitational potential, it is generally 
thought to be unbiased. However, it is a technique that is sen- 
sitive to all the mass along the line of sight, so that projection 
effects may play an important role in adding scatter to any ob- 
served relation. In addition, as it only measures the projected 
(2D) mass, analytical models are needed to transform into the 
more physically motivated spherical (3D) mass, and this is likely 
to add further noise because of cluster triaxiality (e.g., Corless 
& King 2007; Meneghetti et al. 2010). Furthermore, recent the- 
oretical work suggests that some bias may in fact be present 
in weak lensing observations. The systematic 5-10 percent un- 
derestimate of the true mass in the simulations of Becker & 
Kravtsov (201 1) is apparently due to the use of a Navarro-Frenk- 
White (NFW) model that does not describe the data correctly at 
large radii. Notwithstanding, the most recent observational re- 
sults from small samples of clusters for which both X-ray and 
weak lensing data are available indicate either that there is good 
agreement between X-ray and weak lensing masses (Zhang et al. 
2010; Vikhlinin et al. 2009), or that the X-ray mass is system- 
atically lower than the weak lensing mass by up to 20 percent, 
with the underestimate being more important at larger radii (e.g., 
Mahdavi et al. 2008). 

The only investigations of the mass-Fsz relation using 
weak lensing masses published to date have been those of 
Marrone et al. (2009, 2011), using data from the Local Cluster 
Substructure Survey (LoCuSS) 1 . The first directly compared 2D 
quantites (i.e., cylindrical SZ effect vs projected mass) within a 
fixed physical radius of 350 kpc, while the second compared the 
spherically integrated Compton parameter against deprojected 
mass. In the latter case, a much larger scatter than expected was 
found, which the authors attributed to line of sight projection 
effects in weak lensing mass estimates. 

In the present paper we make use of the same weak lensing 
data set from LoCuSS, high quality XMM-Newton archival X- 
ray data, and SZ observations from the Planck 2 All-Sky Survey 
to investigate the interplay between the different mass mea- 
sures and the spherically integrated Compton parameter Ysz in 



1 http : //www . sr . bham . ac . uk/locuss/index . php 

2 Planck (http://www.esa.int/Planck) is a project of the 
European Space Agency (ESA) with instruments provided by two sci- 
entific consortia funded by ESA member states (in particular the lead 
countries France and Italy), with contributions from NASA (USA) and 
telescope reflectors provided by a collaboration between ESA and a sci- 
entific consortium led and funded by Denmark. 



19 clusters of galaxies. We find that for this particular sample, 
the weak lensing mass- Ysz relation at large radii has a slightly 
higher normalisation than that expected from studies based on 
hydrostatic X-ray mass estimates. We show that this is due to the 
hydrostatic X-ray masses being, on average, larger than the cor- 
responding weak lensing masses, in contradiction with the ex- 
pectations from numerical simulations. We show that the prob- 
lem is particularly acute for merging systems and appears to be 
due, at least in part, to a systematic difference in the concen- 
tration as measured by the two methods. In addition, an offset 
between the centres used for the X-ray and weak lensing mass 
determinations appears to introduce a secondary systematic ef- 
fect. 

We adopt a ACDM cosmology with Ho = 70 km s" 1 Mpc" 1 , 

D M = 0.3 and Q A = 0.7. The factor E(z) = V^m(1 + z) 3 + D a 
is the ratio of the Hubble constant at redshift z to its present day 
value. The variables Ma and R& are the total mass and radius cor- 
responding to a total density contrast Ap c (z), where p c (z) is the 
critical density of the Universe at the cluster redshift; thus, e.g., 
^500 = (47r/3)500p c (z)^5 00 . ^ ne Q uan tity Yx is defined as the 
product of Mg 5 5oo, the gas mass within 7?5oo, and 7x, the spec- 
troscopic temperature measured in the [0.15-0.75] Rsoo aper- 
ture. The SZ flux is characterised by Ya, where Y^D 2 h is the 
spherically integrated Compton parameter within and Da is 
the angular-diameter distance to the cluster. All uncertainties are 
given at the 68 percent confidence level. 

2. Sample selection and data sets 

The present investigation requires three fundamental data sets. 
The first is a homogeneously analysed weak lensing data set with 
published NFW mass model parameters to enable calculation 
of the mass at the radius corresponding to any desired density 
contrast. The second is a good quality X-ray observation data 
set that allows detection of the X-ray emission to large radius 
(i.e., at least up to R500). The third is a good quality SZ data 
set including high signal-to-noise SZ flux measurements for all 
systems. 

While there are many weak lensing investigations of individ- 
ual objects in the literature, lensing observations of moderately 
large cluster samples are comparatively rare. For the present 
comparison, we chose to use published results from LoCuSS, 
which is an all-sky X-ray- selected sample of 100 massive galaxy 
clusters at 0.1 < z < 0.3 drawn from the REFLEX (Bohringer 
et al. 2004) and eBCS (Ebeling et al. 2000) catalogues, for which 
gravitational lensing data from the Hubble Space Telescope and 
Subaru Telescope are being accumulated. At the time of writ- 
ing, relevant data from only part of the full sample have been 
published, as detailed below. 

Results from a Subaru weak lensing analysis of 30 LoCuSS 
clusters have been published by Okabe et al. (2010), who provide 
NFW mass model parameters for 26 systems. A similar lensing 
analysis has been undertaken on a further seven merging systems 
by Okabe & Umetsu (2008). Excluding the two of these merging 
clusters that are bimodal and thus not resolved in the Planck 
beam, we have a total of 31 objects for which the best-fitting 
NFW profile mass model is available in the literature. These are 
ideal for our study, given the object selection process (massive 
X-ray clusters) and the fact that the lensing analysis procedure 
is the same for all systems. 

High-quality XMM-Newton X-ray data with at least lOks 
EMOS exposure time are available for 21 of these clusters. Since 
we wished to undertake a fully homogeneous analysis of the X- 
ray data, we excluded two systems, A754 and A2142, whose 
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Table 1: Basic properties of the sample. 
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Notes. X-ray coordinates correspond to the peak of the X-ray emission. The weak lensing coordinates correspond to the position of the BCG. 
Cluster morphological classification is described in Sect. 3.2.4. 



observations consist of a mosaic of several pointings each. For 
the remaining 19 systems, a full hydrostatic X-ray mass analysis 
is possible using the approach described below in Sect. 3.2.2. 



3. Data preparation and analysis 

Table 1 lists basic details of the cluster sample, including name, 
redshift, and the coordinates of the X-ray and weak lensing cen- 
tres. 



SZ observations of the full sample are available from Planck, 
(Tauber et al. 2010; Planck Collaboration I 2011) the third- 
generation space mission to measure the anisotropy of the cos- 
mic microwave background (CMB). Planck observes the sky in 
nine frequency bands covering 30-857 GHz with high sensitiv- 
ity and angular resolution from 31' to 5'. The Low Frequency 
Instrument (LFI; Mandolesi et al. 2010; Bersanelli et al. 2010; 
Mennella et al. 201 1) covers the 30, 44, and 70 GHz bands with 
amplifiers cooled to 20 K. The High Frequency Instrument (HFI; 
Lamarre et al. 2010; Planck HFI Core Team 2011a) covers the 
100, 143, 217, 353, 545, and 857 GHz bands with bolometers 
cooled to 0.1 K. Polarisation is measured in all but the two high- 
est frequency bands (Leahy et al. 2010; Rosset et al. 2010). 
A combination of radiative cooling and three mechanical cool- 
ers produces the temperatures needed for the detectors and op- 
tics (Planck Collaboration II 201 1). Two data processing centres 
(DPCs) check and calibrate the data and make maps of the sky 
(Planck HFI Core Team 2011b; Zacchei et al. 2011). Planck\ 
sensitivity, angular resolution, and frequency coverage make it a 
powerful instrument for Galactic and extragalactic astrophysics 
as well as cosmology. Early astrophysics results, based on data 
taken between 13 August 2009 and 7 June 2010, are given 
in Planck Collaboration VIII-XXVI 2011. Intermediate astro- 
physics results are now being presented in a series of papers 
based on data taken between 13 August 2009 and 27 November 
2010. All of the 19 systems considered in this paper have been 
observed by Planck as part of this survey, and indeed their char- 
acteristics are such that they are almost all strongly detected, 
having a median signal-to-noise ratio of ~ 7. 



3.1. Weak lensing 

As mentioned above, spherical weak lensing masses for the sam- 
ple are given by Okabe & Umetsu (2008) and Okabe et al. 
(2010). These were derived from fitting a projected NFW model 
to a tangential distortion profile centred on the position of the 
brightest cluster galaxy (BCG). In all cases we converted the 
best-fitting NFW profile model to our chosen cosmology and 
obtained Ma by interpolating to the density contrast of interest 3 . 
For the 16 clusters in Okabe et al. (2010), we used the published 
fractional uncertainties at A = 500 and 2500. Uncertainties at 
A = 1000 were obtained from Okabe et al. (2012, priv. commu- 
nication). For the three clusters published by Okabe & Umetsu 
(2008, A520, A1914, and A2034), only M vir and c vir are avail- 
able. Here we estimated the uncertainties at each density contrast 
by multiplying the fractional uncertainty on M v i r given in Okabe 
& Umetsu (2008) by the median fractional uncertainty, relative 
to M v i r , of all other clusters in Okabe et al. (2010) at this density 
contrast. 

3.2. X-ray 

3.2.1. Data analysis 

The preliminary X-ray data analysis follows that described in 
Pratt et al. (2007), Pratt et al. (2010), and Planck Collaboration 
XI (2011). In brief, surface brightness profiles centred on the 



3 For A963, only one-band imaging data are available, which may 
lead to an underestimate of the weak lensing mass (Okabe et al. 2012, 
priv. communication). 
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X-ray emission peak were extracted in the 0.3-2 keV band and 
used to derive the regularised gas density profiles, n e (r), using 
the non-parametric deprojection and PSF-correction method of 
Croston et al. (2006). The projected temperature was measured 
in annuli as described in Pratt et al. (2010). The 3D temper- 
ature profiles, T(r), were calculated by convolving a suitable 
parametric model with a response matrix that takes into account 
projection and PSF effects, projecting this model accounting 
for the bias introduced by fitting isothermal models to multi- 
temperature plasma emission (Mazzotta et al. 2004; Vikhlinin 
2006), and fitting to the projected annular profile. Note that in 
addition to point sources, obvious X-ray sub- structures (corre- 
sponding to, e.g., prominent secondary maxima in the X-ray sur- 
face brightness) were excised before calculating the density and 
temperature profiles discussed above. 

3.2.2. X-ray mass profile 

The X-ray mass was calculated for each cluster as described in 
Democles et al. (2010). Using the gas density n e (r) and temper- 
ature T(r) profiles, and assuming hydrostatic equilibrium, the 
total mass is given by: 



M(< R) 



kT(r) r 
G\im v 



d\wn t {f) d\nT(r) 



dlnr 



dlnr 



(1) 



To suppress noise due to structure in the regularised gas density 
profiles, we fitted them with the parametric model described by 
Vikhlinin et al. (2006) and used the radial derivative d In n Q /d In r 
given by this parametric function fit. The corresponding uncer- 
tainties were given by differentiation of the regularised density 
profile at each point corresponding to the effective radius of the 
deconvolved temperature profile. 

Uncertainties on each X-ray mass point were calculated us- 
ing a Monte Carlo approach based on that of Pratt & Arnaud 
(2003), where a random temperature was generated at each ra- 
dius at which the temperature profile is measured, and a cubic 
spline used to compute the derivative. We only kept random pro- 
files that were physical, meaning that the mass profile must in- 
crease monotonically with radius and the randomised tempera- 
ture profiles must be convectively stable, assuming the standard 
Schwarzschild criterion in the abscence of strong heat conduc- 
tivity, i.e., din T/dlnn e < 2/3. The number of rejected profiles 
varied on a cluster-by-cluster basis, with morphologically dis- 
turbed clusters generally requiring more discards. The final mass 
profiles were built from a minimum of 100 and a maximum of 
1000 Monte Carlo realisations. The mass at each density con- 
trast relative to the critical density of the Universe, Ma, was cal- 
culated via interpolation in the log M - log A plane. The uncer- 
tainty on the resulting mass value was then calculated from the 
region containing 68 percent of the realisations on each side. 

For two clusters, the hydrostatic X-ray mass determinations 
should be treated with some caution. The first is A521, which 
is a well-known merging system. Here the gas density profile at 
large radius declines precipitously, yielding a Jin n e /dln r value 
that results in an integrated mass profile that is practically a pure 
power law at large cluster-centric distances. Although we ex- 
cised the obvious substructure to the north-west before the X- 
ray mass analysis, the complex nature of this system precludes 
a precise X-ray mass analysis. The second cluster for which the 
X-ray mass determination is suspect is A2261, for the more pro- 
saic reason that the X-ray temperature profile is only detected up 
to ^det,max ~ 0.67?5oo ~ 0.87?iooo- In the following, we exclude 
these clusters in cases when the hydrostatic X-ray mass is under 
discussion. 



3.2.3. X-ray pressure profile 

Using the radial density and temperature information, we also 
calculated the X-ray pressure profile P(r) = n Q (r) kT(r). We then 
fitted the pressure profile of each cluster with the generalised 
Navarro, Frank & White (GNFW) model introduced by Nagai 

et al. (2007), viz., 



P(x) = 



feoo, P *) r [ l + feoo, P x) a ] 



(J3-y)/a ' 



(2) 



Here the parameters (a,fi,y) are the intermediate, outer, and 
central slopes, respectively, C5oo, P is a concentration parameter, 
r s = ^50o/c500,p> and x = r/R^oo. In the fitting, the outer slope 
was fixed at ft = 5.49, a choice that is motivated by simula- 
tions since it is essentially unconstrained by the X-ray data (see 
Arnaud et al. 2010, for discussion). The best-fitting X-ray pres- 
sure profile parameters are listed in Table A.l, and the observed 
profiles and best-fitting models are plotted in Fig. A.l. 

3.2.4. Morphological classification 

We divided the 19 clusters into three morphological sub-classes 
based on the scaled central density E(z)~ 2 ft e ,o, which is a good 
proxy for the overall dynamical state (see, e.g., Pratt et al. 2010; 
Arnaud et al. 2010). The scaled central density was obtained 
from a /3-model fit to the inner R < 0.057?5oo region. The 
seven clusters with the highest scaled central density values were 
classed as relaxed 4 ; the six with the lowest values were classed 
as disturbed; the six with intermediate values were classed as 
intermediate (i.e., neither relaxed nor disturbed). Strict applica- 
tion of the R EXCESS morphological classification criteria based 
on scaled central density and centroid shift parameter (w) (Pratt 
et al. 2009) results in a similar classification scheme. 

Images of the cluster sample ordered by E(z)~ 2 ft e ,o are 
shown Appendix B. Henceforth, in all figures dealing with mor- 
phological classification, relaxed systems are plotted in blue, 
unrelaxed in red and intermediate in black. Scaled X-ray pro- 
files resulting from the analysis described below, colour-coded 
by morphological sub-class, are shown in Fig. C.l. 

3.3. SZ 

The SZ signal was extracted from the six High Frequency 
Instrument (HFI) temperature channel maps corresponding to 
the nominal Planck survey (i.e., slightly more than 2.5 full sky 
surveys). We used full resolution maps of HEALPix (Gorski 
et al. 2005) 5 A^de = 2048 and assumed that the beams were 
described by circular Gaussians. We adopted beam FWHM val- 
ues of 9.88, 7.18, 4.87, 4.65, 4.72, and 4.39 arcmin for channel 
frequencies 100, 143, 217, 353, 545, and 857 GHz, respectively. 
Flux extraction was undertaken using the full relativistic treat- 
ment of the SZ spectrum (Itoh et al. 1998), assuming the global 
temperature Tx given in Table 2. Bandpass uncertainties were 
taken into account in the flux measurement. Uncertainties due to 
beam corrections and map calibrations are expected to be small, 
as discussed extensively in Planck Collaboration VIII (2011), 
Planck Collaboration IX (201 1), Planck Collaboration X (201 1), 
Planck Collaboration XI (2011), and Planck Collaboration XII 
(2011). 



4 As these systems have the highest scaled central density, they are 
fully equivalent to a cool core sub-sample. 

5 http://healpix.jpl.nasa.gov 
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Table 2: Masses, X-ray and SZ properties. X-ray masses are calculated as described in Sect. 3.2.2; weak lensing masses were published in the 
LoCuSS weak lensing analysis papers (Okabe & Umetsu 2008; Okabe et al. 2010). 



2500 1000 500 



Name M 250 o D\Y 25 qq M XQm £>i ^iooo M 500 D 2 a Y 500 M gaS;500 T x 7 X ,500 

(10 14 M o ) (10- 5 Mpc 2 ) (10 14 M o ) (10" 5 Mpc 2 ) (10 14 M ) (10" 5 Mpc 2 ) (10 13 M Q ) (keV) (M keV) 



X-ray 



A68 

A209 

A267 

t A291 

t A383 

t A521 

A520 

A963 

A1835 

A1914 

ZwC11454.8+2233 

ZwC11459.4+4240 

A2034 

A2219 

RXJ1720. 1+2638 
$A2261 

RXJ2129.6+0005 

A2390 

A2631 



4 3+0.6 
1.6!°^ 

15^ 
6^ 

i ' U -0.2 

2 9 +o:s 

2> 0+0-4 

5.9!°^ 
4 1+8:^ 

L8 -p| 

9 0+O.6 
-Q-5 
3.3_ 

4.5+°^ 
24+8:§ 

2 5+°^ 

-0.5 

2 5^- 4 
4.9*g 
2.6!i 2 



1:5 



3 ' 8 -Q§ 
2 3 + ° :2 

2.1 

0.9 

5 +m 

1 +01 
39+0.5 

2>1 +o:3 

9 +0 ^ 
5 +0A 

1 9+8:5 

2 9^ 

Z '^-0.3 

7.5 
2.6 
4.9 
1.9 
6.2 



8:§ 



8:5 
-0.4 
3 o +0 - 5 



c Q+0.9 

43^ 

3.0! 05 

2 3 + ^ 
9 9+8:? 
^•^-0.3 
5.6 +2 -° 

6.0 + 



li 
-1.1 

3 9 +0 - 7 
7 7+8:7 
' • ' -0.7 

5 5 +a9 
J -0.8 
9 7+0.5 
7 -0.4 
41+0.8 

41+8:1 

7.2^| 

4 4+0:8 

^ -0.7 

3.1!o°4 1 
3.7t 4 

7.8 + ;- 4 

6.6^ 



1.5 



5.5+JI 

3 6"^ 

8.1!!;I 

8-i-i:8 

40+ - 6 

12^f 
7 1+0:6 

9 0+8:7 

34+0.8 

4 4+0:4 

H -0.4 

14 5 +av 

7.1+8:7 

3.1!°^ 

11 5 + ™ 

7 2 +L1 
'•^-1.1 



6.9!j 2 

6 3 + ' :o 

3.6! - 5 
2 y+8:l 

2 6^ 

^•"-0.4 
an 1+15.0 

+23 
4 9+1.0 

6.9 + ;- 4 

a +0.6 



3.4 
5.4 



+¥:? 
+oi 



10.3 



5.7 

3:8 
6 0+^ 

u.u_ 12 

3 9 +0 - 6 

4 3 4 - ^ 
9.1 + _\l 
9 8 +3 - 8 



:>+0:9 



6.5 + ;i 
10.7 
4.5 
1.9 
0.9 

14.1 + | 
10.6 + ' 

5.3 + f* 
14.7 + 

8.5 + ; 

3.7 + ; 

4.5 
5.5 
21.2 
6.5 
9.0 
3.9 



+°:§ 



15.7 + J 

9.9 + ; 



7 9+0.1 
40+o-o 

o 0+8:8 

-Q.Q 
o 7+O.O 
J - 7 -0.0 

6 5 +ao 

U ' J -0.0 

11 4 +01 
11 -^-o.i 
6>7 +o.o 

1 1 6 + ™ 

10.8+ , 
49+0:0 

7.0^ : ° 

7 + ° :o 
'• u -o.i 

17.6! 02 

6>9 +8:i 

94+8:2 

7 0+8:0 

is.si; 



$8 

u - 7 -0.1 
5.5+0.1 

4 + ^ 
-Q4 
4 9+O.1 

^-0.1 
5-61! 

8-3!;! 

O O+0.2 

8.3_ 02 

4 6 +01 
^• u -o.i 

6 3 +a2 
6 3 +a2 

5-8!:! 

6 7 +a5 

J ' U -0.1 

9 +a2 

^' W -0.2 

7.4^-3 



6>4 +o.3 

7 1 + ^ 

2 2 +ai 

1 5+^ :o 
6+8:8 

o 7+0.I 
J -'-0.1 
Q 1 +0.4 
^• A -0.3 

3 7 +ai 
9 7+0:2 

7 -0.2 

8 9 +a2 

2 3 +ai 

A 4+85 
H -Q.2 

4.4$ 
16.8!- 
4.0!;] 



6.3 + 
4.1! 
14.3! 
7.2; 



3:3 

-0.3 



A68 


1 4+0.6 


A209 


2.1+^ 


A267 


1 4+8:^ 
L - -0.3 


t A291 


9 +04 

W '^-0.4 


t A383 


1 7+0.2 


t A521 


1 1+ 8: 2 


A520 


i ' i -0.4 


A963 


1 o +0 - 2 


A1835 


liit 


A1914 


i ' U -Q.7 


ZwC11454.8+2233 


9+0.4 


ZwC11459.4+4240 


1 8 + ^ : ^ 


A2034 


1 6+°^ 

i ' U -0.9 


A2219 


3 7+O.6 


RXJ1720.1+2638 


1 o+8:S 


JA2261 


c+0.4 


RXJ2129.6+0005 

A2390 

A2631 


3 l +0 ^ 

9 A+M 
^•^-0.4 



2.3 +0 /5 
2.7 
1.8 
0.7 
0.5 
1.2 
2.1; 
1.4 
6.2 
3.1; 
1.5 
2.0: 

2.i + ; 
6.7 + ; 
2.3 + ; 

5 5 +l 

—i 

i.4 + ; 
4.9 + ; 
3.i + : 



-»+8:t 



i 



2.7 + ^ 
5 l +( 

2.4 + ; 
2.3 + ; 
2.6 + ; 
2.4 + ; 

2.5 + | 
2.4 + | 

6.o + ; 

3.1 + 

i.7 + ; 
2.9 + ; 

3.4 + 

6.o + ; 
2.9 + ; 
5.9 + ; 
2.9 + ; 
5 1 +( 

3.7 + ! 



Weak lensing 




4 4+1-0 
^•^-1.0 


4 1 +L2 


7 3 + ® 

a 9+8:§ 


^-j.O 

3 2^ - 7 


15^ 

i ' J -0.7 


4.0+^ 


0.8!;' 


3 3 + ° ) " 7 ' 


44+0.7 


q+8:^ 

D ' y -Q.7 


c 0^8:7 


4.1+11 


34+0.5 
11 2+^ 


4 9+0:9 

^•^-0.7 
9 - 5 ^5 


5.4 +0 - 5 
2 6^ : ^ 


4J+i:6 


2.6 ; ' : ^ 


36^ 
41^ 
13 4+8i 


o+?:8 
-0.9 
5.1!;i 

8.0!|;5 


4 O 4 "^^ 

^' W -0.4 


3.7 +L1 


8> 4+o.8 


8^ 


9 q+8:7 
^"^-0.7 


4 6+i-i 
^• u -i.o 


9>8 +o.7 
6 3 +0 ^ 


7.0 +L3 
4 8 +o:7 



6 +1 3 

U ' W -1.3 

12.3+} 
4.3: 1 
2.1+j 

i.o + £ 
8.0 
9.r 

c 9+6:8 

i4.i + ; 

3.5! 

5.o + ; 

5.4+° 
19.6! 
5 5 +0 

—0 

io.4 + ; 

4.0 + | 

i4.3 + ; 

9.0 + ! 



8.0!^ 

10.6!° 

41+0. 

3.7+0: 
3 7!° 
6.6!| 
11.3! 
6-7!^; 

11.6!; 

10.9!° 

49+0: 

7.0; 

7.0_ 

17.6!* 

7.0! 
92+0 

73!| 
15.8! 
9.8!* 



;8:1 



8.3 + o- 3 

6.6!; 

5.6!| 

3.9 + ; 
4.1!; 
6.r; 
7.9 + ; 
5.6 + ; 
8.4 + ; 
8.5 + ; 
4.6 + ; 
6.4 + ; 
6.4 + ; 
9.6 + ; 
5.9 + ; 
6.r; 
5.6 + ; 
9.i + ; 
7.5 + ; 



6.6!0; 3 
6 9+0.2 

9 0+8:1 

^ -Q.l 

1 4+0.1 
i -^-o.i 
1.5! ao 

4 0!^ : ' 

9.0® 

3.8!S:i 

9 7+0.2 
y - 7 -0.2 
9 2 +0 - 2 
9 0+8:0 



4.5 

4.4 

5 6 +a8 
4i^8:f 

H,A -Q.l 

144+0.5 
7 4^8:4 

7 '^-0.3 



Notes. T x is the spectroscopic temperature within ^500. 

tA291 and A383 were excluded from scaling relation fits involving SZ quantities (see Sect. 3.3 for details). 

$A521 and A2261 were excluded from scaling relation fits involving hydrostatic X-ray mass estimates (see Sect. 3.2.2 for details). 



We extracted the SZ signal using multi-frequency matched 
filters (MMF, Herranz et al. 2002; Melin et al. 2006), which op- 
timally filter and combine maps to estimate the SZ signal. As 
input, the MMF requires information on the instrumental beam, 
the SZ frequency spectrum, and a cluster profile; noise auto- 
and cross-spectra are estimated directly from the data. The al- 
gorithm can be run in a blind mode, where the position, normal- 
isation and extent are all determined by the MMF (e.g., Planck 
Collaboration VIII 201 1), or in a targeted mode, where the posi- 
tion and size are estimated using external data, and only the nor- 



malisation (or SZ flux) is determined by the MMF (e.g., Planck 
Collaboration XI 2011). Here we adopted the latter mode, us- 
ing the position, size, and SZ profile of each cluster determined 
from external X-ray and/or weak lensing data. In this case the 
MMF thus returns only the integrated SZ flux and its associated 
statistical uncertainty. 

Our baseline SZ measurement involved extraction of the SZ 
flux from a position centred on the X-ray emission peak using 
the observed X-ray pressure profile of each cluster described 
above in Sect. 3.2.3 as a spatial template. Apertures were de- 
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RCo[mj NCjny mZ [ivy 

Fig. 1: Relations between 7 S z and total mass for apertures determined from weak lensing mass profiles corresponding to density contrasts of 
A = 2500 (left), 1000 (middle), and 500 (right). In all panels the dark grey region represents the best-fitting relation obtained with slope and 
normalisation as free parameters, and the light grey region denotes the best-fitting relation obtained with the slope fixed to the self-similar value 
of 5/3. Previous results from Marrone et al. (2011), Bonamente et al. (2008), and the analysis of 62 nearby systems by Planck Collaboration XI 
(2011) are shown for comparison. The masses in the latter two studies were derived from X-ray analyses. The original cylindrically integrated SZ 
signal measurement in Bonamente et al. has been converted to a spherically integrated measurement assuming an (Arnaud et al. 2010) profile. 



Table 3: best-fitting parameters for the weak lensing mass-Z)^ 7 50 o 
scaling relations. 



A 


M 
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V Y\M 
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(%) 
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-4.53 + 0.04 


1.48 + 0.21 
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Notes. Relations are expressed as E(z) 7 [D 2 A Y 500 ] = 10 4 [E(zf M/M ] B , 
withy = -2/3, k= 1. 

cr ± is the orthogonal dispersion about the best-fitting relation. 
o~ y\m is the dispersion in Y at given M for the best-fitting relation. 



termined independently either from the weak lensing or the X- 
ray mass analysis. The extraction was achieved by excising a 
10° x 10° patch with pixel size If 72, centred on the X-ray (or 
weak lensing) position, from the six HFI maps, and estimating 
the SZ flux using the MMF. The profiles were truncated at 5 R^m 
to ensure integration of the total SZ signal. The flux and corre- 
sponding error were then scaled to smaller apertures (R500, ^1000, 
^2500) using the profile assumed for extraction. 

We undertook two further tests of the SZ flux extraction 
process. First, we measured the SZ flux using the "univer- 
sal" pressure profile as a spatial template. Here we find that 
the error- weighted mean ratio between these measurements and 
those using the X-ray pressure profile as a spatial template is 
^GNFw/^univ = 1.02 + 0.05, with no trend with morphological 
sub-class. Second, we measured the SZ flux with the position 
left free. In this case the mean error- weighted ratio between the 
flux measurements is Ff ree / ^fix = 1 .04+0.05, again with no trend 
with morphological sub-class. 

The SZ flux measurements for two systems are suspect. One 
object, A291, appears to be strongly contaminated by a radio 
source. The other, A3 83, while not obviously contaminated by 
a radio source and appearing to be a very relaxed system in X- 



rays, exhibits an offset exceeding 4 arcmin (~ 0.87?5oo) between 
the SZ and X-ray positions. This cluster is detected at a rather 
low signal-to-noise ratio, and we also note that Marrone et al. 
(2011), in their Sunyaev-Zeldovich Array (SZA) observations, 
found that this system has an unusually low SZ flux for its ap- 
parent mass. In addition, Zitrin et al. (2012) find that A383 is 
a cluster-cluster lens system, where the nearby z = 0.19 clus- 
ter is lensing a more distant z - 0.9 object that lies about 4' to 
the north-east from the centre of the main system. Zitrin et al. 
also mention the presence of at least two other well-defined op- 
tical structures within 15' of A3 83. In view of the complexity of 
these systems, we exclude them from any analysis that follows 
involving discussion of the SZ signal. 

4. Results 

4.1. Fitting procedure 

We obtained the parameters governing scaling relations between 
various quantities by fitting each set of observables (X, Y) with a 
power law of the form 

E(z) y Y= l0 A [E(z) K (X/X )] B , (3) 

where E(z) is the Hubble constant normalised to its present-day 
value and y and k were fixed to their expected self- similar scal- 
ings with redshift. The fits were undertaken using linear regres- 
sion in the log-log plane, taking the uncertainties in both vari- 
ables into account, and the scatter was computed as described 
in Pratt et al. (2009) and Planck Collaboration XI (2011). The 
fitting procedure used the BCES orthogonal regression method 
(Akritas & Bershady 1996). In addition to fitting with the slope 
and normalisation free, we also investigated the scaling relations 
obtained with the slope fixed to the self- similar values. All un- 
certainties on fitted parameters were estimated using bootstrap 
resampling. 

4.2. SZ-weak lensing mass scaling relation 

Figure 1 shows the relation between the weak lensing mass 
and the SZ flux D 2 A Y^oq measured using our baseline method. 
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Fig. 2: Comparison of present SZ flux measurements to our previous results. Quantities are measured within the R 5 oo derived from weak lensing. 
Left panel: Relation between D 2 A Y 500 and C X sz ix,500 - ^g,500?x, where T x is the spectroscopic temperature in the [0.15 - 0.75] R 50 o region. The 
grey shaded region is the best-fitting power-law relation obtained with slope fixed to 1 ; the red line shows the results from our previous analysis 
of 62 local systems (Planck Collaboration XI 201 1). Right panel: Correlation between D 2 A Y^qq and Mj* derived from the relation of Arnaud et al. 
(2010), compared to the results from Planck Collaboration XI (2011). The shaded region illustrates the best-fitting BCES orthogonal regression 
and associated ±\cr uncertainties. 



All quantities have been integrated in spheres corresponding to 
A = 2500, 1000, and 500, as determined from the weak lensing 
mass profiles. The best-fitting power-law relations are overplot- 
ted, both for regression with the slope fixed to the self-similar 
value of 5/3 (light grey region), and for regression with the slope 
and normalisation free (dark grey region). Numerical values for 
the best-fitting relations, including the dispersion about them, 
are given in Table 3. 

For fits where the slope and normalisation were left as free 
parameters, the slope of the M™ L -D^ Y500 relation is compatible 
with the self-similar value of 5/3 at all values of the density con- 
trast A. The orthogonal scatter about the best-fitting relation is 
cr ± ~ 30-35 percent, and the scatter in D 2 A Y500 for a given M^ L 
is o~ y\m ~ 20 percent, with no significant trend with density con- 
trast. 

Similar fits to the D 2 h Y5oo-M^ L relation for different SZ ex- 
tractions, e.g., with the "universal" pressure profile, or with the 
SZ position left as a free parameter, did not yield results signifi- 
cantly different from those described above. 

5. Discussion 

5.1. Comparison to previous results 

For the D 2 A Y5oo-Mj^ L relation, our results are in good agree- 
ment with earlier determinations at all values of A, albeit within 
the relatively large uncertainties of both our analysis and those 
of previous investigations. A comparison to the most recent re- 
sults of Marrone et al. (201 1), who also use weak lensing masses 
shows that, while the normalisations are in agreement, the slopes 
are slightly (although not significantly) shallower. This is easily 
explained by our exclusion of A3 8 3 from the regression analy- 
sis (see Sect. 3.3 for details); this object was not excluded in the 
Marrone et al. regression fits (see discussion in their Sect. 4.3). 



A fit of their data excluding A383 yields a slope of 1.77 ± 0.16, 
in good agreement with our value (Marrone 2012, priv. commu- 
nication). 

The scatter we observe (cty\m ~ 20 percent) is also in excel- 
lent agreement with that seen by Marrone et al. (201 1). Although 
numerical simulations predict that there is intrinsically only of 
order ten percent scatter between the mass and the integrated 
Compton parameter (e.g., da Silva et al. 2004), observational 
measurement uncertainties and complications due to, e.g., mass 
along the line of sight or cluster triaxiality introduce a further 
source of scatter. Simulations that take an observational ap- 
proach to measurement uncertainties (Becker & Kravtsov 201 1) 
predict a dispersion of order 20 percent, as observed. 

Perhaps the most interesting outcome from the present anal- 
ysis concerns the normalisation of the D\ Ysqq-M^ relation. 
As shown in the right-hand panel of Fig. 1, there is a normal- 
isation offset when the slope is fixed to the self- similar value 
of 5/3, which is significant at 1.4<x with respect to our previ- 
ous investigation of 62 local (z < 0.5) systems using masses 
estimated from the M5oo-Ix,50o relation (Planck Collaboration 
XI 2011). Interestingly, a similar offset was found by Marrone 
et al. (2011) when comparing their scaling relations to those of 
Andersson et al. (2011). This normalisation offset could be due 
either to a larger than expected SZ flux, or to a difference in mass 
measurements between studies (or indeed, both effects may con- 
tribute). Notably, both Planck Collaboration XI and Andersson 
et al. used cluster masses estimated from the Msoo-^x^oo rela- 
tion calibrated using X-ray observations. In the following, we 
first verify the consistency between the present results and our 
previous work, finding excellent agreement. We then examine in 
more detail the reasons behind the observed normalisation offset. 
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Fig. 3: Left panel: Relation between the mass derived from the hydrostatic X-ray analysis (Sect. 3.2.2) compared to the mass derived from the 
^50o-Jx,5oo proxy relation of Arnaud et al. (2010). The shaded region shows the best-fitting BCES orthogonal fit to the data and associated ±\cr 
uncertainties, and the dashed line denotes equality. Right panel: Relation between the mass derived from the hydrostatic X-ray analysis and the 
weak lensing mass of Okabe & Umetsu (2008) and Okabe et al. (2010). The shaded region is the best-fitting regression between the two quantities 
with the slope fixed at 1 . The dashed line denotes equality. 



5.1 .1 . SZ measurements 

Here we wish to verify the consistency between the present 
SZ flux measurements and those we published in Planck 
Collaboration XI (2011). We compare the measurements in a 
statistical sense since not all of the present sample appear in the 
62 ESZ clusters published in that paper. We first compare the 
SZ flux measurement to its X-ray analogue lx,500- This quantity, 
which was first introduced in Kravtsov et al. (2006), is defined as 
the product of the gas mass and the temperature. For consistency 
with our previous work, we define Ix,500 = ^g,500?x, where Tx 
is the spectroscopic temperature in the [0.15 - 0.75] R500 region. 
The left-hand panel of Fig. 2 shows the relation between the SZ 
flux and Ix,500» The latter has been normalised by 

Cxsz = ^L-L = 1.416X10--^- (4) 

m e c z jd e m p M kev 

for yu e = 1.148, the mean molecular weight of electrons for 
a plasma of 0.3 times solar abundance. The grey shaded area 
shows the best-fitting power-law relation between the two quan- 
tities obtained with the slope fixed to l 6 . For comparison, 
we plot the Cxsz Ix,500-^a ^ 500 re l a ti° n obtained by Planck 
Collaboration XI (201 1): D\ F 50 o/C X sz *x,500 = 0.95 ± 0.04. As 
can be seen, the present SZ flux measurements are in excellent 
agreement with our previous determination. 7 

We recall that in Planck Collaboration XI (2011) the mass 
was estimated from the M5oo-^x,500 relation of Arnaud et al. 
(2010). As a second test, we thus calculated Mj* for all ob- 
jects using the Arnaud et al. relation and compared the result- 

6 A fit with the slope left free is shallower, but compatible with unity 
at the 1-sigma level. 

7 The ratio is also consistent with that predicted solely from 
REXCESS X-ray observations: D\ Y 500 /C X sz *x,soo = 0.924 + 0.004 
(Arnaud et al. 2010). 



ing correlation between D 2 A Y^oq-M^ to our previous measure- 
ments. The right-hand panel of Fig. 2 shows this comparison, 
where the grey shaded area is the best-fitting power-law rela- 
tion between the two quantities obtained using orthogonal BCES 
regression with the slope and normalisation as free parame- 
ters. The D\ Ysoo-M^ relation from our previous investigation 
(Planck Collaboration XI 2011) is overplotted. Once again, the 
SZ flux measurements are in very good agreement with our pre- 
vious determination. 

Given the excellent agreement with previous results, we thus 
conclude that the normalisation offset in the D 2 h Y^-M^ rela- 
tion is not due to a systematic overestimation of the SZ flux with 
respect to our previous measurements. 

5.1.2. Mass measurements 

The normalisation offset of the D 2 h Y^-M^ relation may also 
be due to a systematic difference in mass measurements. We first 
need to verify that the hydrostatic X-ray mass estimates detailed 
in Sect. 3.2.2 (M™ ) are in agreement with the expectations from 
the mass proxy relation (Mj* ). The comparison between these 
two quantities is shown in the left-hand panel of Fig. 3. The 
shaded region enclosed by the BCES orthogonal regression fit 
and its uncertainties is entirely consistent with equality between 
the two quantities. 

This leaves us with only one remaining possibility to explain 
the normalisation discrepancy in the D 2 A Y500-M500 relation: a 
systematic difference in X-ray and weak lensing masses. The 
right-hand panel of Fig. 3 shows the comparison between the 
hydrostatic X-ray mass M™ and the weak lensing mass 
A clear offset can indeed be seen. 

However, contrary to expectations, the offset indicates that 
on average the hydrostatic X-ray masses are larger than the 
weak lensing masses. A power-law fit with the slope fixed to 
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Fig. 4: Left panel: This plot shows the ratio of weak lensing mass to hydrostatic X-ray mass as a function of the ratio of NFW mass profile 
concentration parameter from weak lensing and X-ray analyses. Right panel: The ratio of weak lensing mass to hydrostatic X-ray mass is a 
function of offset between X-ray and weak lensing centres. In both panels, the solid line is the best-fitting orthogonal BCES power-law relation 
between the quantities, and the different sub-samples are colour coded. 



1, denoted by the grey region in Fig. 3, indicates that = 
(0.78 ± 0.08) M™. In other words, for this sample, the weak 
lensing masses are ~ 20 percent smaller than the hydrostatic X- 
ray masses at the 2.6<x significance level. 

The mass discrepancy is clearly dependent on morphologi- 
cal sub-class. For relaxed systems, a power-law fit with the slope 
fixed to 1 yields a mean ratio of = (0.94 ± 0.10) M™, 

indicating relatively good agreement between weak lensing and 
X-ray mass estimates. In contrast, the mean ratio for the interme- 
diate and disturbed systems is Af^ = (0.72 ± 0.12) M™. So the 
mass discrepancy is essentially driven by the difference between 
the hydrostatic X-ray and weak lensing masses of the interme- 
diate and disturbed systems (although there is still a slight offset 
even for relaxed systems). 

Zhang et al. (2010) compared hydrostatic X-ray and LoCuSS 
Subaru weak lensing data for 12 clusters, finding excellent 
agreement between the different mass measures [M^ = (1.01 ± 
0.07) M™]. As part of their X-ray-weak lensing study, Zhang 
et al. (2010) analysed the same XMM-Newton data for ten of 
the clusters presented here. For the clusters we have in common, 
the ratio of X-ray masses measured at Rsoo is M™ mg 500/^500 = 
0.83 ± 0.13. This offset is similar to the offset we find between 
the X-ray and weak lensing masses discussed above, as expected 
since Zhang et al. (2010) found ~ ^500- However for re- 
laxed systems (four in total), we find good agreement between 
hydrostatic mass estimates, with a ratio of Af™ ang 500/^500 = 
1.00 ± 0.09. It is not clear where the difference in masses comes 
from, although we note that for some clusters Zhang et al. cen- 
tred their profiles on the weak lensing centre. This point is dis- 
cussed in more detail below. 



5.2. The mass discrepancy 

Our finding that the hydrostatic X-ray masses are larger than the 
weak lensing masses contradicts the results from many recent 



numerical simulations, all of which conclude that the hydrostatic 
assumption underestimates the true mass owing to its neglect of 
pressure support from gas bulk motions (e.g., Nagai et al. 2007; 
Piffaretti & Valdarnini 2008; Meneghetti et al. 2010). If the weak 
lensing mass is indeed unbiased (or less biased) and thus, on 
average, more representative of the true mass, then one would 
expect the weak lensing masses to be larger than the hydrostatic 
X-ray masses. What could be the cause of this unexpected result? 



5.2.1 . Concentration 

To investigate further, we fitted the integrated X-ray mass pro- 
files with an NFW model of the form 



M(< r) = 4np c (z)S c r 5 s 



ln(l + r/r s )- 



1 + r/r s 



(5) 



where p c (z) is the critical density of the universe at redshift z, the 
quantity r s is the scale radius where the logarithmic slope of the 
density profile reaches -2, and S c is a characteristic dimension- 
less density. This model has been shown to be an adequate fit to 
the mass profiles of many morphologically relaxed systems (e.g., 
Pratt & Arnaud 2002; Pointecouteau et al. 2005; Vikhlinin et al. 
2006; Gastaldello et al. 2007). We emphasise that for the present 
investigation we use it only as a convenient fitting formula that 
allows direct comparison with the equivalent weak lensing pa- 
rameterisations. The best- fitting NFW mass model parameters 
are listed in Table D.l and plots of the integrated mass profiles 
and the best-fitting models can be found in Fig. D.l. 

Figure 4 shows the ratio of weak lensing to hydrostatic X-ray 
mass at R500 in terms of the ratio of the concentration parame- 
ter of each NFW mass profile fit. There is a clear trend for the 
mass ratio to depend on the ratio of the concentration param- 
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Fig. 5: Lensing and X-ray mass profiles for A2631 in /i 7( J kpc (left panel) and in terms of R 5 qo (right panel). 



eters. Indeed, a BCES orthogonal power-law fit to the relation 
yields 



M WL i \ -0.25 ±0.07 

500 = io" - 09±0 - 03 I C500 ' WL ] 



IVI 500 



^500,HE 



(6) 



This result is extremely robust to the presence of outliers in the 
relation and to the radial range used to determine the X-ray NFW 
fit. The median ratio of scale radius to /?5oo,he is r s /7?5oo,HE = 
0.35; excluding the three clusters for which r s /7?5oo,HE > 1 yields 
a slope of -0.25+0. 1 1 . The result indicates that the weak lensing 
analysis finds NFW mass profiles that are, on average, more con- 
centrated than the corresponding hydrostatic X-ray NFW mass 
profiles in disturbed systems. As illustrated in Fig. 5, this in turn 
typically explains the trend for the weak lensing masses to be 
lower than the X-ray masses at R500. 

Recent simulations (and some observations) have found that 
the X-ray "hydrostatic mass bias" is radially dependent (e.g., 
Mahdavi et al. 2008; Meneghetti et al. 2010; Zhang et al. 2010; 
Rasia et al. 2012), presumably due to the ICM becoming pro- 
gressively less virialised the further one pushes into the cluster 
outskirts. The difference in concentration that we find here can- 
not be solved by appealing to such a radially dependent X-ray 
"hydrostatic mass bias". If this effect is real, then the hydro- 
static X-ray mass estimates effectively ignore it, meaning that at 
each radius at which the X-ray mass profile is measured, the true 
mass would be underestimated and the underestimation would 
become worse with radius. The resulting hydrostatic X-ray mass 
profile would be over-concentrated relative to the true underly- 
ing mass distribution. Correcting for this effect would reduce 
even further the measured X-ray concentration, exacerbating the 
effect we see here. 



5.2.2. Centre offsets 

In the present work, the X-ray and weak lensing analyses are 
completely independent, extending even to the choice of centre 
for the various profiles under consideration. We recall that Okabe 
& Umetsu (2008) and Okabe et al. (2010) centred their weak 



lensing shear profiles on the position of the BCG. In contrast, 
our hydrostatic X-ray analysis centres each profile on the X-ray 
peak after removal of obvious sub- structures. 8 The fact that at 
7?5oo the mass ratio vs. concentration ratio seems to be driven by 
the intermediate and disturbed systems (see Fig. 4) suggests that 
the different choice of centre could have a bearing on the results. 

We test this in the right-hand panel of Fig. 4, which shows 
the weak lensing to hydrostatic X-ray mass ratio as a function 
of the offset between the BCG position and the X-ray peak. A 
clear trend is visible, in the sense that the larger the offset Rx-wl 
between centres in units of R500, the larger the mass discrepancy. 
Indeed, an orthogonal BCES power-law fit yields 



10 



-0.27±0.07 



WL 



^500 



(7) 



M WL r 71 1 -0.10+0.04 

m 500 

iKZ 500 

Thus at least part of the difference between X-ray and weak 
lensing mass estimates appears to be due to differences in cen- 
tring between the two approaches. Although the trend is visible 
in each morphological sub-sample, the most extreme deviations 
occur in the intermediate and disturbed systems, which all have 
the largest offsets between X-ray and BCG positions. This is 
a well-known characteristic of the observed cluster population 
(e.g., Bildfell et al. 2008; Sanderson et al. 2009; Haarsma et al. 
2010). 

We tested the effect of using a different centre on the X-ray 
mass for two systems. A2631 displays the largest difference in 
mass ratio as a function of concentration parameter (i.e., it is the 
right-most point in the left-hand panel of Fig. 4) and a moder- 
ate X-ray-weak lensing centre offset ~ 0.14 ^500. A520 exhibits 
the largest difference in mass ratio as a function of X-ray-weak 
lensing centre offset (i.e., it is the right-most point in the right- 
hand panel of Fig. 4), with an X-ray-weak lensing centre offset 
of ~ 0.40 /?5oo. For A2631, the choice of centre does not signifi- 
cantly change either the X-ray mass profile or the parameters of 
the NFW model fitted to it, as can be seen in Fig. 5. However, 
when the X-ray profiles were centred on the weak lensing centre 



8 This is in fact required, since otherwise the X-ray analysis would 
give unphysical results. 
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we were unable to find any physical solution to the hydrostatic 
X-ray mass equation (Eq. 1) for A520. 

We note that the dependence of the mass ratio on centre shift 
is qualitatively in agreement with the results of the simulations 
by Rasia et al. (2012). These authors found that the strongest 
weak lensing mass biases (with respect to the true mass) oc- 
curred in clusters with the largest X-ray centroid shift, w. This 
conclusion is supported by the clear correlation between (w) and 
Rx-wl, for which we obtained a Spearman rank coefficient of 
-0.70 and a null hypothesis probability of < 0.001. 

5.2.3. Other effects 

Several other effects could systematically influence the weak 
lensing mass measurements of Okabe & Umetsu (2008) and 
Okabe et al. (2010), and thus contribute to the offset between 
X-ray and weak lensing masses that we find here. 

Firstly, a potential bias could arise from dilution of the mea- 
sured weak lensing signal produced by any cluster galaxies con- 
tained in the galaxy samples used to measure the tangential dis- 
tortion profiles. As discussed in Okabe & Umetsu (2008) and 
Okabe et al. (2010), for all but one cluster considered here 
(A963), data in two passbands were available, enabling a separa- 
tion of cluster and background galaxies based on their location 
in a colour-magnitude diagram. The sample of assumed back- 
ground galaxies consisted of two components: a "red sample" 
with (depending on the available Subaru data) V - V - Iq, 
V - Rq, or g' - Rq colours significantly greater than the colour 
index of the red sequence formed by early-type cluster galax- 
ies; and a "blue sample" with significantly lower colour index 
than the red sequence. The red sample should have very little 
contamination, as all normal galaxies redder than the observed 
red sequence would be predicted to lie at higher redshifts than 
the cluster. By contrast, the blue sample will be contaminated at 
some level by cluster dwarf galaxies undergoing significant star 
formation. Interactions with other galaxies and the intrac luster 
medium in the central regions of the cluster would tend to de- 
stroy dwarf galaxies or quench their star formation, producing 
an observed radial distribution of dwarf galaxies which is much 
shallower than that of the bright early-type galaxies (Pracy et al. 
2004). Hence, while likely present, a small residual contamina- 
tion of cluster galaxies in the "blue" background galaxy sample 
would be difficult to identify and remove without adding more 
photometric filters to the data set. A dilution bias systematically 
lowering the measured weak lensing masses by a few percent 
cannot be excluded, and there would also be significant cluster- 
to-cluster variations in the strength of this effect. 

Secondly, the derived cluster masses are sensitive to the es- 
timated redshift distribution of the gravitationally lensed back- 
ground galaxies. For the lensing mass measurements of Okabe & 
Umetsu (2008) and Okabe et al. (2010), the background galaxy 
redshifts were estimated from the photometric redshift catalogue 
of galaxies in the COSMOS survey field (Ilbert et al. 2009). 
Given the depth of the Subaru data, mass measurements of clus- 
ters in this redshift range (0.1 < z < 0.3) are not very sensi- 
tive to percent level uncertainties in the photometric redshifts 
of the background galaxies. However, a potentially significant 
bias may arise from effectively excluding gravitational lensing 
measurements of galaxies that are smaller than the point spread 
functions (PSFs) of the Subaru images without imposing a sim- 
ilar size cut on the COSMOS galaxy catalogue. This would also 
tend to lower the lensing mass estimate by removing the ap- 
parently smallest (and thus on average more distant) galaxies, 



resulting in an overestimate of the mean effective background 
galaxy redshift compared to the true value. 

Thirdly, the weak lensing distortion measurements of Okabe 
& Umetsu (2008) and Okabe et al. (2010) are based on an imple- 
mentation of the KSB+ method (Kaiser et al. 1995; Luppino & 
Kaiser 1997; Hoekstra et al. 1998). Tests against simulated lens- 
ing data (Heymans et al. 2006; Massey et al. 2007) indicate that 
this method is generally affected by a multiplicative calibration 
bias, underestimating the true lensing signal by up to 15 per- 
cent, depending on details in the implementation of the method. 
This would result in an underestimate of the cluster masses by 
the same amount; however, tests of the particular KSB+ imple- 
mentations of Okabe & Umetsu (2008) and Okabe et al. (2010) 
against realistic simulated weak lensing data would be required 
to measure the calibration factor needed to correct for this effect. 

Finally, recent simulations predict that at large radii the true 
mass distribution departs from the NFW model that was used to 
fit the tangential shear profile in Okabe & Umetsu (2008) and 
Okabe et al. (2010). For example, Oguri & Hamana (2011) re- 
cently investigated the use of smoothly truncated NFW profiles, 
finding them to be a more accurate description of their simu- 
lated clusters. They predicted that the use of a standard NFW 
model would produce an overestimate of the concentration, and 
corresponding underestimation of mass, of order five percent for 
clusters similar to those studied here. This bias can be minimised 
by including the effect of large-scale structure in the shear mea- 
surement uncertainties, thus giving lower statistical weight to the 
weak lensing measurements at large radii (e.g., Hoekstra 2003; 
Dodelson 2004; Hoekstra et al. 201 lb). 

Interestingly, while all these effects are individually well 
within the statistical errors of the mass measurement of a sin- 
gle cluster, they would all have a tendency to bias the measured 
weak lensing masses downwards with respect to their true value. 
Hence, their cumulative effect (combined with the centre offsets) 
could go a long way towards explaining the difference between 
X-ray and weak lensing masses. 

5.3. Summary 

For the present sample, the mass discrepancy between hydro- 
static X-ray and weak lensing mass measurements is such that, 
at 7?5oo, the X-ray masses are larger by ~ 20 percent. This ap- 
pears to be due to a systematic difference in the mass concen- 
tration found by the different approaches, which in turn appears 
to be driven by the clusters that were classified as intermedi- 
ate and/or disturbed. Such systems also show a clear tendency 
to have smaller weak lensing to X-ray mass ratios as a function 
of offset between the X-ray peak and the BCG position. On the 
other hand there is a mass normalisation offset for relaxed sys- 
tems, but it is not significant, as found by previous studies (e.g., 
Vikhlinin et al. 2009; Zhang et al. 2010). 

The results discussed above bring to light three fundamen- 
tal, interconnected obstacles to a proper comparison of X-ray 
and weak lensing mass measurements. The first is how to define 
the "centre" of a cluster and determine what is the "correct" cen- 
tre to use. For X-ray astronomers, the obvious choice, once clear 
sub- structures are excluded, is the X-ray peak or centroid (these 
are not necessarily the same). In contrast, since the weak lens- 
ing signal is most sensitive on large scales, many weak lensing 
analyses use the position of the BCG. 

For relaxed clusters, the BCG position and the X-ray peak 
generally coincide, so that the question of what centre to choose 
does not arise. However, as numerical simulations will testify, 
for disturbed clusters neither of these choices is necessarily the 
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"correct" centre, in the sense that often neither coincides with the 
true cluster centre of mass. Indeed, Hoekstra et al. (201 la) show 
that the recovered weak lensing mass at a given density contrast 
depends on the offset of the BCG from the cluster centroid and 
that the mass can be underestimated by up to ten percent for rea- 
sonable values of centroid offset. This means that for unrelaxed 
systems with offsets between the position of the X-ray peak and 
the BCG position, both approaches will likely give incorrect re- 
sults. 

The second obstacle is connected to the fact that X-rays mea- 
sure close to 3D quantities, while lensing measures 2D quanti- 
ties, and an analytical model is required to transform between 
the two. The choice of an NFW model is motivated by its sim- 
ple functional form and from the fact that it provides a relatively 
good fit to both X-ray and weak lensing data. But we should 
not forget that the original NFW model profile was defined for 
equilibrium haloes (Navarro et al. 1997), and many subsequent 
works have shown that the functional form is not a particularly 
good description of non-equilibrium haloes (e.g., Jing 2000). 
Indeed, very recent work by Becker & Kravtsov (201 1) and Bahe 
et al. (2012) shows that the use of an NFW model profile in a 
weak lensing context can introduce non-negligible biases into 
the mass estimation procedure, primarily because of a departure 
of the mass distribution from the NFW form at large cluster- 
centric radii. 

This point is exacerbated by the third obstacle: the fact that 
X-rays and weak lensing observations have fundamentally dif- 
ferent sensitivities to the mass distribution in a cluster. X-rays 
probe the central regions, and with present instruments at least, 
it is difficult to make very precise measurements of the loga- 
rithmic temperature gradient at and beyond R500. In contrast, 
this is just the radial range at which weak lensing starts to be- 
come most sensitive. Certainly, the combination of strong and 
weak lensing offers much tighter constraints on the inner mass 
distribution (e.g., Kneib et al. 2003), but precise strong lensing 
measurements are difficult to achieve with ground-based instru- 
ments. 

There are several requirements for future progress on this i s- 
sue. The good agreement found here between mass estimates for 
relaxed systems is encouraging and needs to be confirmed for 
a much larger sample of objects, allowing a more precise ob- 
servational constraint to be put on the "hydrostatic mass bias". 
Additionally, such a sample of relaxed objects will allow con- 
straints to be put on the irreducible scatter between the differ- 
ent mass estimates due to the conversion from 2D to 3D quanti- 
ties. Projects such as the "Cluster Lensing and Supernova Survey 
with Hubble" (CLASH; Postman et al. 2012) will surely make 
great progress on these questions. However, a new, or at least 
co-ordinated, approach is needed in the case of dynamically dis- 
turbed systems. Here, numerical simulations can also be used to 
inform the different analyses and to optimise the mass estima- 
tion procedure in each case. One possible approach would be to 
use the centroid of the projected X-ray pressure profile as the 
point around which both X-ray and weak lensing profiles could 
be centred. 

Eventually, X-ray and weak lensing mass estimates will be 
needed for a representative (or complete) sample of systems. We 
stress that data quality across such a sample must be as close to 
homogeneous as possible. In the case of the X-ray data set, the 
data must be sufficiently deep to measure the temperature profile 
to 7?5oo. Similarly stringent data quality will also be required for 
the weak lensing data set. 



6. Conclusions 

A well-calibrated relation between the direct observable and 
the underlying total mass is essential to leverage the statistical 
power of any cluster survey. In this paper we presented an in- 
vestigation of the relations between the SZ flux and the mass 
for a small sample of 19 clusters for which weak lensing mass 
measurements are available in the literature and high-quality 
X-ray observations are available in the XMM-Newton archive. 
This "holistic" approach allowed us to investigate the interde- 
pendence of the different quantities and to attempt to square the 
circle regarding the different mass estimation methods. 

Using weak lensing masses from the LoCuSS sample (Okabe 
& Umetsu 2008; Okabe et al. 2010), we found that the SZ flux 
is well correlated with the total mass, with a slope that is com- 
patible with self- similar, and a dispersion about the best- fitting 
relation that is in agreement with both previous observational 
determinations (Marrone et al. 2011) and simulations that take 
into account observational measurement uncertainty (Becker & 
Kravtsov 201 1). However, at R500, there was a normalisation off- 
set with respect to that expected from previous measurements 
based on hydrostatic X-ray mass estimates. 

We verified that the SZ flux measurements and hydrostatic 
mass estimates of the present sample are in excellent agreement 
with our previous work (Planck Collaboration XI 2011). The 
normalisation offset is due to a systematic difference between 
hydrostatic X-ray and weak lensing masses, such that for this 
particular sample the weak lensing masses are 22 + 8 percent 
smaller than the hydrostatic X-ray mass estimates. The differ- 
ence is essentially driven by the intermediate and morpholog- 
ically disturbed systems, for relaxed objects, the weak lensing 
mass measurements are in good agreement with the hydrostatic 
X-ray estimates. 

We examined the possible causes of the mass discrepancy. 
At 7?5oo, the X-ray-weak lensing mass ratio is strongly corre- 
lated with the offset in X-ray and weak lensing centres. It is 
also strongly correlated with the ratio of NFW concentration 
parameters, indicating that the mass profiles determined from 
weak lensing are systematically more concentrated than the cor- 
responding X-ray mass profiles in disturbed systems. We argued 
that a radially dependent "hydrostatic mass bias" in the X-ray 
observations would exacerbate this effect, and discussed several 
other alternative explanations, including dilution and uncertain- 
ties due to the use of NFW mass profiles to model the weak 
lensing data set. 

Significant progress on the mass calibration of clusters can 
only be achieved with a dedicated X-ray-lensing survey of a rep- 
resentative sample of clusters. Data of sufficient quality are an 
essential prerequisite for such a survey. X-ray observations with 
sufficiently deep exposures to measure the temperature at 7?5oo 
are needed, as are optical observations of uniformly high quality. 
We expect such data to become available in the coming years. 
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Appendix A: Pressure profiles and best-fitting 
model 

The X-ray pressure profile of each cluster was fitted with a 
generalised-NFW (GNFW) model as described in Sect. 3.2.3. 
Table A. 1 gives the best- fitting model and Figure A. 1 shows each 
X-ray profile with the best-fitting GNFW model overplotted. 

Appendix B: Image gallery 

Figure B.l shows the 0.3-2keV band X-ray image gallery of 
the cluster sample, arranged from top left to bottom right in 
order of the morphological characterisation parameter n Q $, the 
central density. Images are corrected for surface brightness dim- 
ming with z, divided by the emissivity in the energy band, taking 
into account galactic absorption and instrument response, and 
scaled according to the self-similar model. The colour table is 
the same for all clusters, so that the images would be identical if 
clusters obeyed strict self- similarity, and each panel is 1.257?5oo 
on a side. 

Appendix C: Scaled X-ray profiles 

Figure C.l shows the X-ray profiles of the sample. They have 
been radially scaled by the Rsoo determined from the M500- 
^x,5oo relation of Arnaud et al. (2010, see also Pratt et al. 2010). 
Relaxed (or equivalently, cool core) systems are plotted in blue, 
disturbed systems in red, and intermediate objects in black. The 
gas density profiles are scaled by the expected self- similar evo- 
lution with redshift; the temperature profiles are scaled by the 
average spectroscopic temperature in the [0.15 - 0.75] R500 re- 
gion; the pressure profiles are scaled by P500; the mass profiles 
are scaled by M500 (also estimated from the M500 - ix,soo rela- 
tion). 

Appendix D: Mass profile fits 

The integrated mass profile of each cluster was fitted with an 
NFW model as described in Sect. 5.2. Table D.l gives the cor- 
responding best- fitting NFW model and Fig D. 1 shows the mass 
profile of each cluster with the best-fitting NFW model overplot- 
ted. 
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Table A.l: best-fitting X-ray pressure profile parameters. 
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Notes. Column (2): R500, the radius corresponding to a density contrast of 500, estimated from the weak lensing mass analysis of Okabe & 
Umetsu (2008); Okabe et al. (2010). Column (3): P 500 as defined by Eq. 5 of Arnaud et al. (2010). Columns (4) to (7) give the best-fitting GNFW 
parameters for the pressure profiles (Eq. 2). The external slope parameter ft has been fixed to 5.49 (see text). No uncertainties are given as the 
GNFW parameters are highly degenerate. The profiles and best-fitting GNFW models are illustrated in Fig. A.l. 



Table D.l: best-fitting mass profile model parameters. 
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Notes. Columns (3,5): concentration parameter. Columns (4,6) M 50 o 
from the best-fitting NFW model. The profiles and best-fitting X-ray 
NFW models are illustrated in Fig. D.l. 
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Fig. A.l: Pressure profiles of the sample with the best-fitting GNFW model overplotted (red line). The dotted vertical line indicates R 5 oo for each 
cluster. 
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Fig. B.l: Image gallery. 
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Fig. C.l: X-ray profiles scaled by R 5 oo from the M 50 o-ix,5oo relation of Arnaud et al. (2010). From left to right, top to bottom: gas density, gas 
temperature, gas pressure, integrated mass. Cool core systems are plotted in blue, morphologically disturbed objects in red, and intermediate 
systems in black. The grey line in the pressure plot is the universal pressure profile of Arnaud et al. (2010). 
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Fig. D. 1 : X-ray mass profiles of the sample with the best-fitting NFW model overplotted (red line). The dotted vertical line indicates R 5 oo for each 
cluster, determined from the best-fitting NFW model. 
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